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impulse approximation. The reaction allows for the study of the production 
process in the nuclear medium without being obscured by the details of nu- 
Cji I clear transition densities. First, a pion electroproduction operator suitable for 

nuclear calculations is obtained by extending the Blomqvist-Laget photopro- 
duction operator to the virtual photon case. The operator is gauge invariant, 
unitary, reference frame independent, and describes the existing data reason- 
ably well. Then it is applied in nuclei to predict nuclear cross sections under a 
I variety of kinematic arrangements. Issues such as the effects of gauge-fixing, 

the interference of the A resonance with the background, sensitivities to the 
quadrupole component of the A excitation and to the electromagnetic form 
factors, the role of final-state interactions, are studied in detail. Methods on 
how to experimentally separate the various pieces in the coincidence cross 
section are suggested. Finally, the model is compared to a recent SLAC ex- 
periment. 
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I. INTRODUCTION 



The goal of this work is to develop a theoretical framework for analyzing the exclu- 
sive, quasifree pion electroproduction from complex nuclei, denoted by A{e, e'7TN)B, in the 
A(1232) resonance region. Here, "exclusive" means that the outgoing particles e', tt, N 
are detected in coincidence, and the nucleus undergoes transitions to discrete final nuclear 
states, usually the ground state. The term "quasifree" refers to those processes which can 
be identified as taking place on a single nucleon inside the nucleus (impulse approximation). 

The study of this reaction is of interest because it provides a testing ground for our 
knowledge on several areas of nuclear physics research: the electromagnetic production of 
pions from nucleons; quasi-elastic electron scattering off nuclei; the pion-nucleus and the 
nucleon-nucleus interaction. The combination is best appreciated by visualizing the photo- 
production process in the impulse approximation with a single-particle model of the nucleus. 
In this picture (see Fig. |I]), the incident photon (real or virtual) penetrates the nucleus and 
couples to an individual nucleon via the latter 's charge and magnetic moment. This causes 
the nucleon to oscillate and radiate pions. The produced pions, along with the nucleons 
that also exit the system, subsequently rescatter from the remaining nucleons before finally 
escaping and reaching the detector. The initial pion production stage involves knowledge of 
the photoproduction of pions off single nucleons, the rescattering stage requires an under- 
standing of the way pions and nucleons scatter off nuclei, while the overall transition of the 
nucleus between specific initial and final states is the subject of electron scattering studies 
below pion threshold. In principle, each of these ingredients can become the subject of the- 
oretical scrutiny by using the best knowledge on the other ingredients. The most important 
motivation, however, is to use this reaction to study the A excitation and propagation in 
the nuclear medium. 

Traditionally, pion photoproduction from nuclei has been studied with the reaction 
A{'y,7T)B Since this reaction requires the final nucleon to remain bound in the resid- 
ual nucleus, it involves relatively high momentum transfers. Consequently, these processes 
are very sensitive to the details of the nuclear transition densities which often obscure the 
process of primary interest: pion photoproduction in the medium. In recent years, it is rec- 
ognized that sensitivities to the nuclear transition matrix elements can be greatly reduced 
by allowing the final nucleon involved in the pion production process to exit the nucleus, 
namely, by studying the coincidence reaction ^(7, 'kN)B. This is mainly due to the quasifree 
nature of the reaction: the momentum transfer to the residual nucleus can be made small. 
By measuring in coincidence the decay products of the A, one can get a much better handle 
on the A in the medium. Comparison with existing data already showed good promise of 
such reactions 0-0. 

Pion electroproduction from nuclei in the A region, on the other hand, has been studied in 
less detail. While the reaction A(e, e'n)B has been investigated in Ref. 0, to our knowledge 
there is no theoretical work on the process A{e, e'7rN)B in the A region. The use of electrons 
instead of photons adds the advantage of probing the longitudinal response and longitudinal- 
transverse interference, but with the price of dealing with more complicated structure in the 
cross section. 
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Such coincidence experiments are challenging because the final nuclear state, the outgoing 
pion and nucleon all need to be identified in coincidence with sufficient energy resolution 
and solid angle coverage. Data for A{'y,TiN)B in the A region are sparse 0. Recently, two 
experiments on the reaction have just been completed and data are being analyzed [^H- 
Data on A{e, e'iTN)B are almost non-existent, until the first experiment was recently carried 
out at SLAC . However, with the advent of new high duty cycle accelerators such as LEGS 
at Brookhaven, Bates, NIKHEF, MAMI and CEBAF, the situation is expected to greatly 
improve in the near future. 

In section ||, we describe our pion electroproduction operator and show how it compares 
with existing data. In section |IT|, the DWIA formalism for the reaction A{e, e'7TN)B is 
derived. In section |^ we show results under a variety of kinematic arrangements in order 
to expose different aspects of the physics involved, and compare with the SLAC experiment. 
Section gives a summary of our findings. The full operator is given in Appendix The 
electromagnetic form factors used in the calculation are given in Appendix 



II. PION ELECTROPRODUCTION ON THE NUCLEON 



In order to study pion electroproduction in the resonance region on complex nuclei, one 
first has to understand the production process on the nucleon. To this end there have 
been extensive studies, for example, in |[lO|-pO|. The expressions for the elementary cross 
sections have become standard. We have rederived them as a by product in our derivation 
of the DWIA formalism for the nuclear case pO| using the density matrix method. This also 



serves as a consistency check. In the following, we will skip the derivation and only give the 
expressions relevant for our discussions. 



A. Differential cross section 



All quantities are in the laboratory frame unless otherwise mentioned (cm. quantities are 
explicitly denoted by a superscript c). Due to the relative weakness of the electromagnetic 
interaction, electron scattering can be treated as the exchange of a virtual photon which 
carries energy uj = — Ee' and momentum k = pe — Pe'. The 4-momentum transfer 
squared is spacelike and given by 

e = uj^-\k\'^ = -AE^Ee' sin^ — (1) 

where 6e is the electron scattering angle and the electron mass is neglected. The differential 
cross section for unpolarized pion electroproduction on single nucleon can be written as an 
electron ffux factor times the virtual photon cross section : 



dVl(,idEf^idVlT^ di^TY 
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where 

dfl-rr dQjr ' dQjT ■ \j \ ■ / 
The flux factor is defined by 



day dar , da^ , darr , /——, — Z da^^ 



27r2 Ee 1 - e ^ ^ 

In this expression, a is the fine structure constant, K = {W'^ — m'^)/2m^ where m is the 
nucleon mass, is called the virtual photon equivalent energy and W"^ = k'^ + rn^ + 2muj is the 
invariant mass squared of the vrN pair. The variable e is the degree of transverse polarization 
of the virtual photon and is given by 

/ 2|kP 2 0eV^ 

Note that e is invariant with respect to the Lorentz boost along the virtual photon direction 
k. The transverse, longitudinal, polarization and interference cross sections, as they are 
called, are given in terms of the hadronic current matrix elements by 

^ = 1 E (iwp+iwr), (6) 



da,, , —k"^ 



E mw (7) 



d^TT 2 m^m^ 



^ = ^V-^ E iiJ'^riJ^) + w(J^Y) (9) 

where the z-axis is along k. Here the notation (J^) = (m/| \mi) represents the matrix 
elements between initial and final nucleon states. Eq. (|]) to Eq. @ are valid both in the 
laboratory frame and in the cm. frame. The kinematic factor A is given in the laboratory 
frame by 



m q 




?>2t^'^K\En + 


(q- k) ■ q/|q| 
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^ = oo,.^lz. , J7L „.„I2P (10) 

and in the cm. frame by 
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Ik 






2K 





Ik^l V47rVr; ^ ^ 



where the pion four momentum = (i?^,q). The above definitions ensure that in the 
real photon hmit /c^ — *• 0, the transverse cross section Eq. (H) goes to the real photon cross 
section. 

Two comments are in order. First, a different form is sometimes used in the literature 
where the explicit pion angular dependence is written out: 

Inr^.^'^.^'ln:'''' ^.cos(20.) + Ve(l + e)^sm^.cos0.. (12) 

Under this definition the polarization and interference cross sections will differ from those in 
Eq. @ by the extra angular factors. Therefore one has to be careful about the definition used 
when comparing with data or with other calculations. It is obvious from this expression that 
the total virtual cross section has contributions only from the transverse and longitudinal 
terms: 

o"^ = o"t + ecTi- (13) 

Second, we discuss a kinematic singularity problem in the cm. frame. The virtual 
photon energy in the cm. frame can be written as uj'^ = {W^ + k"^ — m'^)/{2W). Since 
the virtual photon is spacelike, u'^ can vanish or even become negative. This would lead to 
apparent singularities in the cross sections (see Eqs. (|^) and (|])). For pion electroproduction, 
W > m + ttLt^ ~ 1080 MeV, so it could happen for 4-momentum transfer starting at —k"^ = 

+ 2mm^ ~ 0.24 {GeV/cY which lies mostly in the region of interest. In the A region, 
for example, W=1230 MeV, then uj^ vanishes at approximately —k"^ = 0.63 [GeV/c)'^. One 
obvious way to circumvent the problem is to work in the laboratory frame. However, since 
all the electroproduction data are conventionally presented in the cm. frame, we still need 
to transform the cross sections into the cm. frame. This can be achieved by the following 
Jacobian: 

:i7^ = 1 - ^/?cos^^ (14) 

dn% |q| |q| 

where the Lorentz boost from the laboratory frame to the cm. frame is /3 = k/(c(j + m). 
We will give another solution below when we discuss gauge invariance. 



B. The electroproduction operator 

Traditionally, models of pion electroproduction on the nucleon are built in terms of 
CGLN amplitudes |[T0|-[12| by means of a multipole decomposition. These amplitudes are 



expressed in the cm. frame. A transformation is needed in order to use them in nuclear 
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calculations where the struck nucleon has a momentum distribution (Fermi motion). Such 
a transformation is complicated and sometimes ambiguous. Moreover, it is not clear how to 
extrapolate the amplitudes to off-shell. 

Our goal is to have an operator that can describe the elementary process reasonably well 
and that is suitable for application to nuclear calculations. To this end, we employed the 
techniques of Blomqvist and Laget (BL) pl| , p2| to construct an electroproduction operator. 



The BL photoproduction operator is based on an effective Lagrangian approach, incorpo- 
rating gauge invariance and unitarity. It is expressed in an arbitrary frame of reference 
which makes it convenient for use in nuclear calculations. The diagrammatic structure of 
the amplitude provides a physically transparent picture of the individual processes which 
contribute and facilitates the discussion of nonlocality and off-shell effects. The BL ampli- 
tudes describe the elementary pion photoproduction data reasonably well over a wide range 
of energies and have enjoyed success in many nuclear calculations These features make it 
a good candidate for our purposes. In recent years, attempts to upgrade the BL model have 
been made |^ and models based on the Hamiltonian approach have been developed |^ . 



The procedure of extending the BL pion photoproduction operator to electroproduction is 
straightforward. Firstly, the same Feynman diagrams as in photoproduction are considered, 
but electromagnetic form factors describing the charge and magnetic moment distribution of 
the particles need to be introduced at appropriate vertices to replace the static charge and 
magnetic moments. They are the pion form factor Fj^lk'^), the nucleon form factors -Fi(/c^) 
and F2{k'^) (or equivalently Guik'^) and GE{k^)), the nucleon axial form factor Fa(A;^), the 
A resonance form factor F^{k'^), and in the case of u meson exchange, the u form factor 
F^{k^). Secondly, terms proportional to e ■ k that were dropped in photoproduction need to 
be included since they are non-vanishing due to the additional longitudinal polarization of 
the virtual photon. Thirdly, in addition to the transverse quadrupole E2 term, there is the 
longitudinal quadrupole C2 term (sometimes called L2). To improve the readability of the 
paper, the full expressions for the obtained operator are given in Appendix ^ and the form 
factors used in the model are given in Appendix 0. 



C. Gauge invariance 



When electromagnetic form factors of the various particles are introduced at each vertex, 
the Born terms in the resulting interaction are no longer gauge invariant (the delta term is 
constructed to be separately gauge invariant and hence not affected). Current conservation 

Jo = |k| Jjuj (15) 

has been assumed when deriving the formalism in the previous section. As the magnitude 
of fc^ increases the form factors become quite different from each other, especially the pion 
form factor. Thus, restoration of gauge invariance is needed when different form factors are 
used in the calculation. While various solutions to this problem have been explored, we use 
a method of restoring gauge invariance which does not change the interaction in the Lorentz 
gauge. 
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Let us assume we have written the interaction in the Lorentz gauge as Sf^J^. We can 
check gauge invariance (or more precisely current conservation) by replacing by the photon 
four momentum A;^ and checking if k^J'^ = 0. If contains general form factors then this 
equation will not be satisfied. However, current conservation can be restored by defining a 
new current 



(16) 



Furthermore, in the Lorentz gauge e^J = e^^J^ since the Lorentz gauge condition is e^/c^ = 
0. Restoring current conservation by this method remains ad hoc, but it seems to do minimal 
damage, and allows us to use realistic form factors. With a conserved current we have 
regained gauge invariance and thus can carry out the calculation in the most convenient 
gauge. Once better quality experimental data is available, the validity of this approach can 
be investigated. For the present we will begin with Laget's non-relativistic parameterization 
of pion photoproduction with gauge fixing to the same order in {p/m) and check if our 
different form factors and additional current terms continue to fit the available data. The 
resulting gauge fixing terms in Eq. (0) are given in Appendix 0. 

As discussed above, the vanishing virtual photon energy in the cm. frame is problematic 
since the transformation to eliminate the time component of the hadronic current by current 
conservation becomes singular at this point. It is not a true singularity in the sense that 
one can in principle factorize out u in to explicitly cancel the u factor in front. But in 
practice that is cumbersome and sometimes not possible. 

An alternative is to keep all four components of the current. As a result [Q, the 



transverse cross section darr/dVL^^ and polarization cross section da^r/d^T^ remain unchanged, 
while the longitudinal and interference cross sections take on new forms: 



doj, 
dVL,. 



rriimj 



-k^ 



and 

dcTTL 
dVL^ 



2 ^ 

^ miirif L 



{{J.Y{J.) + {J.Y{Jz)) 



\\^\uJ 



{{J.Y{J,) + {JoY{J.)) 



{{JoYW + iJ.YiJo)) 



(17) 



(18) 



Clearly, there is no singularity in this formulation, no matter what reference frame is used. 
Since no current conservation is assumed in the derivation, one can work in any gauge 
that is convenient. One can also use the new equations to investigate the degree of gauge 
invariance violation before gauge fixing. On the other hand, if one uses current conservation 
to eliminate Jo, Eqs. ([T7|) and (|18D will reduce to Eq. (|^) and Eq. (|^), respectively, as they 
should. 



D. Comparison with data 



In this section we compare theoretical calculations on the nucleon performed with the 
electroproduction operator with a large body of existing data p5|-pm, most of which were 



7 



obtained in the 1970s at Bonn and Saclay. Since our aim is to use the operator in nuclei, 
we only select some representative figures to show the overall quality of the operator. More 
details can be found in Ref. [EOl- 



Fig. ^ shows the comparison with the total cross section data on p{e, e') in the A region. 
This experiment was inclusive — only electrons were detected. Thus, the curves represent 
the combined contributions from both the ep — > n7c~^ and the ep — pn^ processes. Since 
the total virtual cross section only consists of longitudinal and transverse contributions (see 
Eq. (p2D), a separation of them is possible by varying electron kinematics. Several interesting 
features can be seen in the figure. First, the A resonance peaks are clearly present in the 
W dependence. Second, the longitudinal cross sections are consistently small compared to 
the transverse cross section. Third, in the real photon limit A;^ — > 0, the calculations agree 
with the measured real photon points. Fourth, the operator is supposed to be valid in 
any reference frame, which is indeed the case judging by the small differences between the 
solid lines and the dotted lines. The overall agreement with the total cross section data is 
satisfactory. 

Fig. shows the separated longitudinal and transverse cross sections with pions exiting 
along the direction of the virtual photon momentum. In such kinematics the polarization 
and the interference cross sections vanish (see Eq. (|12|)), allowing the longitudinal-transverse 



separation by varying e. The solid curves represent the calculation with the full operator 
which uses realistic form factors, contains the E2 and L2 terms and includes the gauge fixing 
terms. The dotted curves are calculated with the same form factors = -Fa = -Fi, in which 
case the Born terms are already gauge invariant, thus there is no need to restore gauge 
invariance. The dashed curves are the same as for the full operator, except the quadrupole 
terms E2 and L2 have been omitted. The dash-dotted curves show the effect of leaving out 
the gauge fixing terms. Clearly, gauge restoration which only affects the longitudinal term 
is important in improving the agreement between the model and the data. 

The remaining figures in this section compare three calculations: the full operator (solid 
lines), no gauge fixing terms (dotted lines), no E2 and L2 terms (dashed line), in order to 
study the significance of gauge fixing and the quadrupole excitation of the A. The definition 
in Eq. (p!^ ) is used for the polarization cross section daTr/d^n and the interference cross 
section daTL/dfl-,, in these figures. 

Fig. ^ and Fig. |, compares the W dependence of the separated cross sections with 
data for different values of k"^ and 6^^. Gauge fixing is not very significant in these cases, 
although its inclusion improves the agreement between theory and experiment. There is 
some sensitivity to the E2 and L2 terms around W=1230 MeV, right on top of the A 
resonance. 

In Fig. ^, the pion azimuthal angle distributions of the virtual cross sections (unsepa- 
rated) are presented along with data at different values of W, k'^, and 9-„. Note that both 
gauge fixing and the quadrupole excitation exhibit more significance in the A region. 

Fig. 1^ shows the result of an analysis of the p{e, e'n^)p reaction performed at Bonn, in a 
nearly coplanar (^^r = 20*^) and a perpendicular {(p^^ = 90°) kinematic condition. The strong 
difference between these two kinematics is due to the polarization cross section of which the 
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contribution changes sign when 0^ varies from to 90''. In perpendicular kinematics, the 
contribution of the polarization cross section vanishes and the experiment determines the 
transverse coupling of the A. In nearly coplanar kinematics the effects of the polarization 
cross section are at maximum and the quadrupole component becomes more important. 
Note, that the effects of gauge fixing in this case are small, almost indistinguishable from 
the full calculation. 

Fig. § presents the separated cross sections at the even higher value of fc^ = —1 [GeV/cY. 
The agreement with data in this case is excellent. 

The overall agreement agreement between the model and the data is satisfactory, which 
is sufficient for our purposes in this exploratory study. We stress that we have not attempted 
to readjust any parameters in order to fit the data. All the parameters are already fixed by 
the BL pion photoproduction operator. In this sense the results are predictions rather than 
fits. As new data become available in the near future, it may become necessary to refine 
the operator by refitting, or to construct a new operator suitable for use with a relativistic 
treatment of the nucleons in the nucleus. 



III. THE DWIA MODEL FOR PION ELECTROPRODUCTION ON NUCLEI 

Having obtained the elementary transition operator, we are now in the position to 
study quasifree pion electroproduction on nuclei, A{e, e'TTN)B, in the A region. The term 
'quasifree' means that the process can be identified as taking place on a single nucleon in 
nuclei. This happens when the missing momentum (or the momentum transfer to the target) 
is relatively small, say, less than 300 MeV/ c. In such a situation, the impulse approximation 
is expected to hold. We employ the Distorted Wave Impulse Approximation (DWIA) which 
closely follows our previous approach for photoproduction in Ref. 0. However, the new lon- 
gitudinal polarization of the virtual photon generates additional terms in the cross section 
which are nontrivial in a full DWIA framework. See Fig. |l] for an diagrammatic illustration 
of the reaction A{e, e'TTN)B in the A region in DWIA. We first discuss the kinematics, then 
briefly describe the nuclear structure inputs, and finally derive the cross sections. 

A. Kinematics 

In this section, we define our coordinate system and discuss some kinematic aspects of the 
exclusive reaction (e, e'nN) on nuclei. In the laboratory frame, the four-momentum of the 
incoming electron is Pe = {E^, Pe), of the outgoing electron p^' = {E^', Pe'), thus the virtual 
photon carries the four-momentum k = Pe — Pe' = (i^^k). The outgoing pion and nucleon 
have four-momenta of q'^ = (i?^, q) and p'^ = {E^, p), respectively. The target nucleus is at 
rest with mass Mi and the recoiling residual nucleus of mass Mf has the momentum 

Q = k-q-p (19) 
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and kinetic energy Tq = The momentum transfer Q is also referred to as the missing 
momentum. Overall energy conservation requires 

u; + M, = E^ + EN + Mf + TQ. (20) 

As shown in Fig. the z-axis is defined by the virtual photon direction k, and we choose 
the azimuthal angle of the pion, 0,^ = 0, by defining y = k x q and x = y x z. This puts the 
X-axis in the electron scattering plane while the y-axis is normal to the electron scattering 
plane. In general, the pion and the nucleon can both go out of the plane. 

We assume that the reaction takes place on a single bound nucleon with four-momentum 
p^i = {Ei, Pi) and that energy and momentum are conserved at this vertex [i.e., the impulse 
approximation). Thus pi = — Q and Ei = E.„ + E^ — oj. If Q does not vanish, the struck 
nucleon is off its mass shell. This is the only reasonable off-shell choice since the photon, 
the pion, the outgoing nucleon, and the recoiling nucleus are all external lines and must be 
on their respective mass shells. 

The magnitude of the momentum transfer to the recoiling nucleus has a wide range, 
including zero, depending on the directions of the outgoing pion and nucleon with respect 
to the virtual photon. However, since the reaction amplitude is proportional to the Fourier 
transform of the bound state single particle wave function, it becomes quite small for mo- 
mentum transfers greater than about 300 MeV/c. Thus for all but the lightest nuclei we 
can safely neglect the nuclear recoil velocity (and recoil energy Tq) and generate optical 
potentials for the outgoing particles in the laboratory frame. As noted in |Q , since the resid- 
ual nucleus can take up varying amount of momentum but little recoil energy, the reaction 
offers great kinematic flexibility and by appropriate choices one can investigate the produc- 
tion operator, the bound state wave-function, or the final state interaction of the outgoing 
meson and nucleon. 

Another useful kinematic quantity is the invariant mass of the outgoing vrA^ pair W = 
\J{p + q)"^ which indicates whether the pion production process takes place in the A(1232) 
resonance region. Note that there is a difference in kinematics between the real and virtual 
photon cases. In our case, the magnitude |k| is always larger than u, while in the real case 
they are equal. This results in different off-shell behaviors of the production amplitudes. 

B. Nuclear Structure Inputs 

As discussed in the previous section, we are primarily interested in cases of low momen- 
tum transfer to the recoiling nucleus. Thus, our choice of single particle wave-functions for 
the bound state is not critical as long as the basic size of the orbital is described correctly. 
For convenience we use harmonic oscillator wave-functions which have the advantage that 
their Fourier transforms are simply obtained. For each nucleus under consideration, we 
adjust the harmonic oscillator range parameter b until the RMS radius of the ground state 
charge distribution agrees with the experimentally determined values. 
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For the continuum nucleon wave-functions we solve the Schrodinger equation with an 
optical potential present whenever the particle is to be detected. Further, we use the ex- 
perimentally determined separation energies for a given orbital in order to fix the value of 
the mass of the recoiling nucleus Mf. Many optical models for the outgoing nucleon are 
available. We use a non-relativistic reduction of the global optical model of Clark and col- 



laborators This model has the advantage that it fits nucleon scattering over a wide 

range of energy and A values, and hence is very useful for making surveys of a wide range of 
possible experimental situations. Other models have been tried |^ and were found to make 



little difference. Once experimental data is available for the exclusive reaction, an optical 
model specific to the nucleus and energy range of the outgoing nucleon can be substituted. 



For the continuum pion wave-functions, the SMC optical model [32| is used. Alterna- 



tively, we have tried the pion global optical potential in which is based on multiple 
scattering and found that the two approaches give similar results. 

Another nuclear structure input is the overlap of the residual nucleus with the A-1 
spectator particles which just provides an overall normalization, or spectroscopic factor, to 
the cross section. These have been taken from quasi-elastic electron scattering studies. 



C. Differential Cross Sections 



The differential cross section for y4(e, e'TiN)B in the laboratory frame can be written as 

, Ef. nie nie dpe' da uIm dv} Mf dO 

d(T — — 

IPel E,,{2Tif 2E^{2nf E^ {271^ Ef{2Tif 

X Y.\Mn? (27r)45(4) (p, + -p,, - - g -p) (21) 

where denotes the sum over final spins and the average over initial spins and the (5-function 
defines the overall energy- momentum conservation. After integrating over Q and |p| and 
summing over the nuclear part of the transition matrix element, the 5- fold coincidence 
differential cross section can be written as a sum over single particle matrix elements times 
some kinematic phase space factor: 



d^a 


|Pe 


\ml 


MfrriN 




|P 




dEe' dfle' dEj, dVL^, dVL tv 


Pel 


(27r)3 2{2Tif\EN + Ef - ENV- 


(k - q)/p2 



X 2 ^ YTi\T{s.s\m,,a)\\ (22) 

s s' 
nis a 

Here, a = {nljm} is the quantum number of the bound nucleon, s and s' are electron spins, 
rris is the spin projection of the outgoing nucleon and Sa is the spectroscopic factor. The 
single particle transition matrix element is given by 

T(s,s',m„«) = J rf3^vl>(+)(r,-p)0(+)(r,-q)j^^J'^e^'^-'-vl/,(r) (23) 
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In this equation, = u{pe', s') 'jf^u^pe, s) is the electron current, e the electron charge, 
"^a the bound nucleon wave-function, "^^^ , 0^"^-* the distorted outgoing nucleon and pion 
wave-functions, and J'^ the hadronic transition current. Introducing an auxiliary current 

fl/. = jfi -—k^ (24) 



and using the fact k^J'^ = and = 0, one can eliminate the time component of the 
hadronic current and work with the space components only, i.e., jfj^J'^ = —aiJi = —ai6ijJj = 
— J2x(^^^iW ^jWJj where we have inserted a dyadic to project the expression onto the 
virtual photon polarization basis which is chosen to be unit vectors along the three coordinate 
axes (A = x, y, z). Under this basis, the matrix element squared with sums over the electron 
spins can be written as 

Em' = |^Ep-'<^v (25) 

ss' XX' 



where we have defined the electron density matrix elements 

P.v = E • a]* [e{X') ■ a] (26) 

ss' 



and the hadronic matrix elements 

w,ims,a) = J rfV vI/W(r, -p) ^^^r, -q) s(A) ■ J e^^'^ ^^(r). (27) 



The matrix element Wx{ms,a) has the identical form to that in the real photon case (see 
Eq. (6) of Ref. Q]) except here one has the longitudinal polarization in addition to the 
two transverse polarizations. Therefore, the method for evaluating Eq. (^) closely follows 
Ref. and we will not be repeated here. 

After carrying out the electron spin sums, and with the help of Eq. (|I|) and the relations 



Pe 



Pe' 



|Pe 




Pe'l 




k 





Sin^e, P^=Pe 



0, 



IPe 



|Pe/| COS^^e 



Pe 



[-\Pe'\ + |Pe| COS 



(28) 



we find the 3x3 electron density matrix 



-k' 



1 



P 



ml 1 — e 



1 /-2fc2 







^e(l + e) 




#^(l + e) 







uj2 t / 



(29) 



where e has been defined in Eq. (|^). Using this matrix, the 5-fold differential cross section 
Eq. (p2D can be cast into the form of an electron flux factor times the virtual photon cross 
section, just as in the single nucleon case: 
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dEe' d^l^i dEj^ d^l-jj^ d^l jv dE-,^ dil-,^ d^ljy 

except that the electron flux factor here is defined by 

r = J^E^ ^7 1 /o-i \ 

27r2 Ee 1 - e 

where ii^^ = (s^ — )/2Mj, and s = /c + Pj. The virtual photon cross section is given by 

d^a„ 



dE-n dQ.,r dQ 



N 



d^ar + ed'^ai^ + ed'^aTT + Je{l + e) d^arL, (32) 



where we have introduced the short-hand notations for the transverse cross section: 

"'--raa^-^S^^sO-i^-KD. (33) 

the longitudinal cross section: 

- dE^ dn^ dn^ - ^ 2j + 1 co^ ^^''^ 

the polarization cross section: 

^ ^ s, 



d' 



a. 



dE^ dn^ dn^ 2j + 



cT. ^i(Kr-Kr), (35) 



and the interference cross section: 



,3 _ i^aTL „ ^ S„ -1 /-2fc2 



The kinematic factor C is given by 

^ _ MfUiN |q| IpI , . 



The above definitions of the factors F and C ensure that in the real photon limit —k'^ —>■ 0, 
the transverse cross section d^aT reduces to the real photon cross section (see Eq. (3) of 
Ref. i). 
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IV. RESULTS 



In this section we present our results for exclusive quasifree pion electroproduction on nu- 
clei. There exist almost no data for the reaction, except for one experiment recently carried 
out at SLAC which will be discussed below. However, the experiment ^^0{e,e'7i~py^0 
is being planned at the new high duty cycle electron accelerator at Mainz |^ . Such exper- 
iments are difficult to perform since they require large solid angle acceptance, good energy 
resolution, and long beam time. We expect our results to be helpful in planning such ex- 
periments. We use ^^C as the generic target and consider processes leaving the residual 
nucleus in its ground state. This corresponds to knocking out a p-state nucleon from ^^C. 
Our formalism is equipped to study other nuclei as well. Our goal is to get an overview 
of the magnitude and shape of the cross sections and to seek out suitable kinematics for 
future experiments. For this particular purpose, it is sufficient to use the plane wave impulse 
approximation (PWIA) which costs much less computer time. The role of distortions will 
be discussed separately. 

We present our results using three kinematic arrangements: 

I) the magnitude of the missing momentum is kept fixed (quasifree kinematics), 

II) the magnitudes of all momenta are fixed (fixed kinematics), 

III) the missing momentum is allowed to vary freely (open kinematics). 

We will consider coplanar setups first, i.e., the pion plane and the nucleon plane coincide 
with the electron scattering plane with = 0° and = 180° (see Fig. |). Out of plane 
arrangements will be presented when we discuss the separation of cross sections. 



A. Results for Kinematics I 



In this kinematic arrangement, we keep u, k, and the magnitude of the missing momen- 
tum Q fixed, and study the pion angular distribution. We refer to it as 'quasifree kinematics' 
since it most resembles the free two-body kinematics. It is achieved in two steps. For a given 
pion angle, we set Q = and calculate the corresponding 'two-body' kinematics inside the 
nucleus. Note that this is not exactly the same as the free two-body case since nuclear 
binding has to be taken into account. Next, we select a finite value for Q (or Pm), thus 
violating momentum conservation (Eq. ([T9|)). However, energy conservation (Eq. (^)) is 
almost maintained because the nuclear recoil energy Tq is small. As a result, the pion and 
proton energies remain at their two-body values for a fixed pion angle, but the proton angle 
changes to a new value. This setup is suitable for studying effects of final state interactions 
since the pion and proton energies both change as the pion angle is varied over the whole 
angular range. 



To be more specific, we show in Fig. |T0| an example for the reaction C(e, eV p) C, 



a- 



which corresponds to knocking out a neutron from the P3/2 orbital. The solved kinematic 
variables (the proton angle 6^ at Q=0, and 6p at finite Q, the proton kinetic energy Tp, the 
pion kinetic energy T^r) are plotted as a function of the pion angle. The invariant mass W of 
the ttN pair remains roughly constant around 1210 MeV over the range. Note that different 
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Q values are necessary to generate relatively large cross sections for orbitals with different 
angular momentum (around zero for s-waves, about 100 MeV/c for p-waves, and about 150 
MeV/c for d- waves, etc). 

Fig. |ll| shows the charged pion virtual cross section along with the individual structure 
functions for such a kinematic setup. The energy of the virtual photon is chosen to be 
uj = AOOMeV, which lies in the A region. The values for u and — fc^ can be varied separately 
by tuning the initial and final electron scattering energies and angles. One can obtain the 
full 5-fold differential cross section, Eq. (pO[), by multiplying the virtual cross section with 
the electron flux factor which depends on the values of electron scattering energies and angle 
actually used in the experiment. In forward pion direction the cross section is dominated 
by the longitudinal structure function. This is similar to the free pion electroproduction 
process on the nucleon where the dominance of the longitudinal term at low t (corresponding 
to small Ott) has been used to extract the pion electromagnetic form factor in a Chew-Low 
extrapolation. Our findings confirm that this situation also occurs in the nuclear case, 
possibly permitting the extraction of the pion form factor in the nuclear medium. At larger 
pion angles, d^a^ rapidly decreases and d^aj. becomes dominant. The interference and 
polarization cross sections contribute moderately for intermediate pion angles with negative 
sign. Shown separately are the contributions from the background Born terms and the 
resonant A term. As they add coherently there are large interference effects between the 
two contributions. Note that the longitudinal structure function is completely dominated by 
the Born terms, more specifically, the main contribution comes from the pion pole term. The 
A term is small for d^aTL whereas in d^a^T the Born terms are small. Thus each structure 
function can provide complementary information on the production process in the nuclear 
medium. 

The quadrupole component of the A has been an important issue in electroproduction 
processes because it relates to the internal structure of the A resonance. Unlike the real 
photon which can only excite the transverse quadrupole component E2, the virtual photon 
allows for the additional longitudinal component L2 (or C2). The quadrupole component is 
small (only a few percent) compared to the dominant magnetic dipole transition (Ml) and 
its value is not well determined. Fig. |l^ shows the sensitivities of the response functions to 
E2 and L2 in neutral pion electroproduction from nuclei by varying the E2/M1 ratio (see 
Appendix We find large sensitivities in the longitudinal and interference terms while the 
transverse cross sections are completely insensitive. This is again similar to the situation on 
the free nucleon. Thus, separating these cross sections would allow to verify if the E2/M1 
ratio is modified in the nuclear medium. These sensitivities appear mostly in the 7r° channel, 
due to the small size of the Born terms. 



Fig. |T3| shows the dependence on the four- momentum squared, k'^, and demonstrates the 
sensitivity to form factors. The solid curve displays the full operator with realistic form 
factors and gauge fixed, the short-dashed curve contains a dipole pion form factor instead 
of a monopole type, the long-dashed curve is the result with all form factors equal. Most 
sensitivities appear in the longitudinal and the interference terms. Note, that under the 
given kinematic conditions k"^ is constrained to below -0.35 (GeV/c)'^. 



Fig. presents the distortion effects due to the final state interactions of the pion and 
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proton with the residual nucleus. For details on the distorted waves employed, see Ref. 
and references therein. One notices the strong energy dependence of the pion distortion 
since the pion energies large for small pion angles and vice versa. The nucleon distortion 
slightly reduces the transverse cross section as well as the longitudinal one at forward angles. 
The DWIA results given in this figure can serve as a guide on how distortion affects the 
PWIA results for kinematics II and III below since they are given at fixed pion and proton 
energies. 



B. Results for Kinematics II 



In this kinematic arrangement, we keep uj, k, the pion energy T^^, the nucleon energy 
T/v and the magnitude Q fixed, and again vary the pion angle. We refer to it as the 'fixed 
kinematics'. Only the angles of the momentum vectors vary. This has the advantage that the 
uncertainties from nuclear structure effects and final-state interactions are minimized, thus 
the production process in the medium can be exposed. Since the lengths of the momentum 
vectors involved are all fixed, the angular range accessible is limited. 



Fig. |T5| shows the contributions from Born and A terms to the charged pion virtual cross 
section and the individual response functions as a function of pion angle. Under the given 
kinematics, the pion angle covers the range from 25° to 120*^. Fig. ^ shows the sensitivities 
to the quadruple moment E2 and L2. 



C. Results for Kinematics III 



For this kinematic situation, we keep u, k, the pion energy T^, and nucleon angle 9p 
fixed. It is equivalent to mapping out the missing momentum pm distribution, we refer to it 
as the 'open kinematics'. Since pm is allowed to vary, the whole kinematic region is opened 
up, although for large Q the cross sections fall off dramatically due to the bound nucleon 
wave-functions. 

Fig. |l^ shows the charged pion angular distributions for one set of such kinematics. The 
twin peaks in the distributions reflect the missing momentum distribution of the struck 
nucleon. In this case it is the P3/2 state in ^^C. There is strong interference between 
the background and the resonance resulting in a forward-backward asymmetry of the two 



contributions in the virtual cross section. Fig. |T8| shows the sensitivities to the quadrupole 
component. Again, large sensitivities are displayed in the longitudinal and the interference 
terms. 

Under these kinematic conditions, when the outgoing pion and proton are both detected 
along the direction of the virtual photon, the polarization and the interference terms van- 
ish. As a result, the longitudinal and transverse pieces can be extracted by a Rosenbluth 
separation. Fig. ITO shows the four- momentum — fc^ distribution and sensitivities to the form 



factors. Large sensitivities to changes in the form factors are displayed in the longitudinal 
term. Note that in this special case, the longitudinal cross section is almost as large as the 
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transverse one. The height of the two peaks are reversed due to the falloff with the form 
factors. 



D. Separation of differential cross sections 

The virtual cross section contains contributions from the transverse and longitudinal 
currents and their interferences, each having different sensitivities to particular informa- 
tion in the production process. It is desirable that these individual terms be disentangled 
experimentally. One may recall that in electroproduction on the nucleon, explicit angular 
dependences can be factorized out (see Eq. (0)), making a separation possible. It would be 
desirable to have a similar separation in pion electroproduction from nuclei. 



In a series of studies on quasi-elastic electron scattering off nuclei [0, Donnelly found 
that angular factors of the azimuthal angles can indeed be pulled out from the response 
functions in the so-called 'super Rosenbluth formula'. Defining the difference 

= 0^ - (t)N (38) 

and the average 

+ 0^) , (39) 



1 
2 



It was found that angular factors of the average can be explicitly separated out, leaving 
the remainder only dependent on the difference A0. Comparing our formalism to that in 
Ref. shows that the same separation scheme is applicable here. Translating into our 



formalism, the angular factors only appear in the polarization and the interference cross 
sections: 

d^aTT = d^cTTT cos 20 + d^a^,^ sin 20, (40) 

d^cFTL = d^a^^ COS0 + d^a^^ sin 0. (41) 

where each term is divided into two terms (with superscripts a and b) that depend only on 
the difference A0 multiplied by the appropriate angular factors that only depend on the 
average 0. 

Mapping out a distribution while keeping A0 fixed allows separating {d^a^ + ed^a^^), 
d^a^T and d^aj.^. Then, one can use the electron scattering angle 6e dependence in e to 
perform a standard Rosenbluth separation of d^aT and d^cTi^. Thus, all cross sections in 
Eq. ( p2|) are in principle experimentally accessible. 

Fig. |20| shows the out-of-plane distribution at A0 = 135° using kinematics III. Also shown 



are the sensitivities to the quadruple moment. Note that the polarization cross section d^aTi 
is very close to a pure cos 20 distribution, indicating that d^a^j. is small compared to d^a^j. 
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Furthermore, the is no dependence in the longitudinal and transverse structure functions. 
The interference term closely resembles a sin function, thus the d^a^^ term is small. We 
found that cross sections drop considerably for configurations of A0 < 90° due to increasing 
missing momentum. Thus, A0 > 90° configurations are kinematically favored. In particular, 
at A(j) = 180°, the polarization cross section is purely cos 20 and the interference cross section 
is purely sin0, as shown in Fig|2l]. At A0 = 0°, d^a-rj. is purely cos 20, and d^axL is purely 
COS0, but the cross sections are several orders of magnitude smaller compared to those at 
A0 = 180°. 



E. Comparison to SLAC data on A{e, e'pir ) 



In a recent SLAC experiment first measurements were made for the reactions (e, e'n~p) 
and (e, e'pp) in the A region on ^H, CO, Ar and Xe targets. The photon-nucleon invariant 
mass W ranged from 1.1 to 2 GeV, and —k"^ from 0.1 to about 1 GeV/c^, with an average 
— fc^ ^ 0.35 GeV/c^. However, due to low statistics, the exclusive differential cross section 
was integrated over the energy and solid angle of the outgoing proton and over the azimuthal 
angle of the outgoing pion. In order to compare with the data, we need to carry out the 
following integration numerically 



rfV r ( d^ 



dEf.1 d^le' d cos On J \ dEe' dVL^.' dE^^ dVt.,^ dVt 



N 



dEndVtNdct)^. (42) 



This integration can be done in PWIA, but becomes immensely time consuming in term of 
computer power for the full DWIA. So far, we approximate the DWIA result by scaling the 
PWIA result with a distortion factor evaluated at free production kinematics (which means 
the missing momentum is zero): 

/ d^a \ fd^anv^TA\ ( d^a \ 

X ,c ,...n (43) 



^ dEe' dVte' d cos 9^ ) ^^j^ \ d^ap^iA J q^q \ dE^' dile' d cos 9^ j p^^^ 



Fig. ^ depicts the pion angular distribution of the integrated cross section in the lab- 
oratory frame relative to the direction of the virtual photon momentum for C0(e,e'p7r~) 
at W=1.210 GeV, and -k'^ = 0.35GeV/c . In obtaining the cross sections for the CO tar- 
get, we have summed over s-shell and p-shell knockout for ^^C and ^^O. The PWIA result 
clearly overestimates the data. The DWIA results roughly agree with the data at backward 
angles, but overpredict at forward pion angles. Fig. ^ also shows that the fall-off of the 
cross section with increasing pion angle is more gradual than predicted in either calculation. 
At forward angles the data are suppressed relative to DWIA by approximately a factor of 
two. This is reminiscent of the results of the Bates experiment |0] for real photons, where 
our calculations along with other studies found good agreement with the data at backward 
angles and were approximately three times larger than the data at forward angles. The 
difference between forward and backward pion angles may be due to the structure of the 
electroproduction operator, which is predominantly the resonant A term at forward angles 
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and non-resonant Born terms at backward angles. This would suggest that the discrepancy 
at forward angles is due to the A-nucleus interaction. However, full DWIA calculations have 
to be carried out before any definite conclusions can be drawn. 

V. SUMMARY 

In this study we have established a DWIA formalism for calculating quasifree pion elec- 
troproduction from nuclei, A{e, e'7TN)B, in the A region. The reaction provides the same 
promise in studying A excitations in nuclei as does the photoproduction reaction ^(7, ttN)B, 
with the added advantage and complexity generated by the longitudinal polarization of the 
virtual photon. The sensitivity to the nuclear structure of the target is minimal. The only 
information required of the target is the single particle bound wave-function, the spectro- 
scopic factor, and the optical potentials. Kinematically, the reaction provides a great deal 
of flexibility since the target can take up a wide range of momentum transfer but little recoil 
energy. 

We obtained the pion electroproduction operator by extending the BL pion photopro- 
duction operator which has enjoyed much success in pion photoproduction from nuclei. It is 

gauge invariant, frame independent, and simple to implement in a non-relativistic descrip- 
tion of nuclei. We compared the operator to a large body of data on the single proton and 
found overall good agreement. 

Nuclear cross sections are predicted under a variety of kinematic situations to help expose 
the different aspects of the reaction. Such issues as the interference between the A resonance 
and the Born background, sensitivities to the quadrupole component of the A excitation, and 
sensitivities to the electromagnetic form factors, are studied. Reasonably large sensitivities 
to the quadrupole component are found in the neutral pion channel. Methods on how to 
separate the various structure functions in the cross section are suggested. These results are 
expected to be useful in planning future experiments. 

The first measurements of the reaction in the A region, a SLAC experiment on 
C0{e,e'p7r^), are analyzed using our framework. Preliminary results show that it re- 
veals a similar pion forward-backward anomaly first found in the Bates experiment on 
^^0{j,pn^y^O. Clearly, more theoretical and experimental studies are needed at this point. 
The SLAC experiment suffers from low statistics and poor energy resolution. As a result, 
the cross sections had to be integrated over, thus making comparisons with calculations 
difficult and ambiguous. However, the situation is expected to improve as more exclusive 
data become available in the near future with machines at MAMI, Bates, LEGS, NIKHEF, 
and the commissioning of CEBAF. 
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APPENDIX A: PION ELECTROPRODUCTION OPERATOR 



In this appendix, we give the full operator for both charged and neutral pion electropro- 
duction. The operator is a straightforward extension of the Blomqvist-Laget pion photo- 
production operator with appropriate form factors and gauge correction terms introduced. 
It is given by t^^Tr = ■ J where is the virtual photon polarization vector and J is the 
pion electroproduction current. We decompose the operator into spin and spin 1 terms by 
writing 



t^,^{X, k, Pi, q, p) = L + i <r • K. 



(Al) 



The non-spin flip term L and the spin flip term K each consists of a coherent sum of the 
Born and A resonance terms: 



L = L 



Born 



(A2) 



K — l^Born + Ka. 



(A3) 



The Born terms in PV coupling for various production channels are given by the following. 

For 7^ -|- p — >• TT"*" -|- n: 



■'Born 



\/2ego 
2m 



q- (k X £a), 



(A4) 



K 



Born 



ego 



2m 



-FAik')+Ffik' 



mE^, 



Ea{PS + Ea) 



ex 



+ 



{q — ky — ml 
G^^ik') GUe) 



2E,{PS - E,) 2E,(P0 - E,) 



2E„(P0 - E,) 
q X (k X £a) > ■ 



(A5) 



For 7^ -|- n 



TT 



■'Born 



p: 

V2ego 
2m 



G%[{k'^ 



+ 



2Ea{PS - Ea) 2E,{P^ - E,) 



q- (k X ea), 



(A6) 



K 



V2 



Born 



ego 



2m 



mE^ 



Eb{P^ + E,) 



ex 



{q-kf 



mt 



+ 



G%ik' 



2E,{PS - E,) 2E,{P^ - E,) 



2Efe(po - E,) 
q X (k X £a) L 



(A7) 
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For 7„ + p ^ 7r° + p: 



T - ^ 

^Born - 



+ 



2m [2Ea{P^ - Ea) 2E,{P^ - E,) 



q 



q) -(kxe;,), (A8) 



'^Born 



2m [[2Ea{PS-Ea) 2E^{P^-E,) 



(q - ^(2Pi - q)) X (k X £x) 



2 £a • (2pi + k) 

'''^'''2E^{P^^-E^) 
^ £x ■ (2p - k) 



E^ 



2Eb(po - E, 
m ' 



m 



Et, 



(2p.-q)-q\ fFf^k'^ + ^pP^k'^ 



2m 



2m 



(A9) 



For 7„ + n ^ 7r° + n: 



■'Born — 



ego 
2m 



+ 



G^{kp- 



2Ea{PS - Ea) 2E,{P^ - E, 



q-|^(2p.-q))-(kx£;,), (AlO) 



K 



Born 



ego 
2m 

+ex 



G^{k'^ 



2E,{PS-Ea) 2E,{P^-E, 



m 



Ea{PS + Ea) 

m 

Eb{Pb° + E,) 



E^ 



Et, 



(2p + q)-q \ f E 

2m 

(2p,^q):q\ (ElLF^n^e^ 
2m / \2m ^ ^ ' 



q-^(2pi-q)) x(kxe;,) 



(All) 



The photon, incoming nucleon, pion and outgoing nucleon four-momenta arc k = (u;,k), 
Pi = {Ep^, Pi), q = {Et,, q), p = (-Eat, p ), respectively. They can be in any reference frame. 
The nucleon mass is denoted as m. The four-momenta in the s and u channels are Pa = k+Pi 
and Pb = Pi — q = p — k and Ea^b = (|Pa,6p + rn^Y^'^. For the n — N coupling constant we 
use go/4:n = 14. 

In the 7r° channels, the following 6i;-exchange term should be added coherently to L. 

1 eg^, 



f^n {q^ — k^^Y — ml 



q • (k X £x) 



(A12) 



where m^ = 750 MeV, g^^ = 10. 

The A resonance terms with Ml and E2 and C2 transitions are given by 



24 



P^-MI + iTaMa 



(A13) 



-eC^a^GiG^FAik'^) f 1 



P^-MI + iTaMa 

(q* ■ k*) ex + (q* ■ ex) k* + 
2exk 



|-q*x(k*X£;,)e^<^'^^ 



2 



m 



[q* ■ k*)k* 



3 m 
2ujma 



l^A ■ Pi) q 



2£Ak 

' 3|k|2 



3(q* ■ k*)k* - |k|'q 



(3Ma + m)(MA + m) 
2 * T 2c(jma 



(3MA + m)(MA + m) |k*| 



(AM) 



where k* = k — (Ma — m)pi/m and q* = q — -E'7r(k + pj)/MA, the isospin coefficients 
C,r±C*7 = and (7,^0(7^ = |. The couphng constants Gi,Gs, the mass of the deha Ma 
and the width of the deha Fa were treated as parameters and fitted to the data. We use 



the following parameterization [21 



Ma 
Ta 



1225 MeV, 



110 



|q| V^^a 1 + (0.007|p„|)^ 



Gi = 0.34 



,\Pa\J Pa l + (0.007|q|)2 

Ma + m 



MeV, 



G3 = 2.18/m^ MeV"^ 



(A15) 



where pa and q are in MeV/c and m,r is in MeV. The phases used to restore unitarity are 
functions of the variable x = Pa — 1080 in MeV and are given in degrees as 



^j^j = -0.1228v^ + 0.0735X, 

i)E = 3.9136x/x + 0.2795a; - 0.00049a;^ 



(A16) 



The symbol a in Eq. ( |A14| ) is a constant that measures the relative strength between the 
Ml and E2 transition amplitudes in the A and here takes the value a = 0.8, while fc'^ = 
(Pf — m^) / {2Pa) is the equivalent photon energy in the vr-N cm. frame. 

In the following we give the gauge corrections for the Born terms. They should be added 
to the Born part of the ex ■ J for various channels. The A term is made separately gauge 
invariant and no additional terms need to be added. 
For 7„ + p ^ 7r+ + n: 



ex 



.jiG) 



Sx 



■k V2 



ego 



-k^ 2m 



I (7 



+ 



{Fa-F^) + ^{FI-Fa)+FI 
2m 



{F. - Fl) - 



^{FI-Fa) q+-(Ff-FA)p. 
2m \ m 

mE^ 



2EaiPS-Ea) EaiP^ + Ea)^ 



(A17) 
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For 7„ + n ^ TT + p: 



j(c)_ ^A-k v/2e^o 



00 



q--{F[-FA)p^ 

m 



+ 



uj ( ijjE mE 



(A18) 



For 7^ + p ^ 7r° + p: 



J(G) 



— /c^ 2m 

2m^-(|kp + 2p,-k)(l + |^) -2m^+(-|kp + 2p.k)(l + ||) 
2Ea{PS-Ea) ^ ' 2E,{P^-E,) 



-Ff 



ojEt, 



ujEtt 



2E,{PS-Ea) 2E,{P^-E, 



+ 
+ 



pEt, 2mu! — |kp — 2pj ■ k j^^p-^Tr —2muj — |kp + 2p ■ k 



Fl 



m 2Ea{PS - Ea) 
E^2mu- |kp-2pi -k 

2Ea{PS-Ea) 

m 



Ea{PS + Ea) 

m 



MPS + E,)\'' 



F? 



m 2E,{P^ - Ek) 

( IjJE^ IjJE. 



^2Ea{PS + E^) 2E,{P^ + E,) 



Err — 



Et, 



(2p + q)-q \ (pv _ E^pP 

2m )\^ 2m ^ 
(2p, - q) • q \ / p E^^p 
2m \ ^ 2m ' 



k . 



(A19) 



For 7„ + n ^ 7r° + n: 



Ett 



E^ — 



2Ea{PS + Ea) V " 

(2p, - q) • q 



2E,{P^ + E,) \ " 2m 



Ett — 



(2p + q)-q ^ 
2m 



(A20) 



Note that these additional terms go to zero (to order (p/m)'^) if the same form factors are 
used: Fa = Ff = F„. The following relations can help demonstrate it: 

2Ea(P° - Ea) ^ 2mu - |kp - 2pi ■ k, 
2Eb(P° - Eb) ^ -2muj - |kp + 2p • k, 
2E,(P° + ^a) ~2m(2m + u;), 

2E,,{PS + Eb) ^ 2m(2m - u). (A21) 



APPENDIX B: ELECTROMAGNETIC FORM FACTORS USED IN THE MODEL 



For the nucleon form factors we use the well-known dipole form: 
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.2\ -2 



1 + 



0.71 



(Bl) 



where /ip = 2.79, /i„ = —1.91, r = — A;^/(4m^) and k"^ is in (GeV/c)^. The form factors Fi 
and F2 are given in terms of and Gm by 



^1 TTT ■ 



1 + r 



(B2) 



(B3) 



Note that in this definition (fc^) = 0. For the pion form factor we use the monopole form 



F.{k') = 1 + 



For the axial form factor we use 



FAik') = 1 + 



0.45 



0.90 



(B4) 



(B5) 



For the A form factor we use the following form 



6/1 0.72 



(B6) 



Finally, the form factor at the ujtt^j vertex is taken to be of p(775) type: 



F^{k^) = 0.374 1 + 



0.60, 



(B7) 



The RMS charge radius •\/< > is related to the form factor by (assume A;^ < 0): 

dF{e] 



< -6 



d{-k^) 



(B8) 



fe2=0 



27 



FIGURES 




FIG. 1. An illustration of the reaction A(e,e'nN)B in the A region. 
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FIG. 2. The fc^ and W dependence of the separated total cross sections for the inclusive reaction 
p{e, e') is compared to the data [^]. The upper curves are for the transverse cross section ax, lower 
curves for longitudinal cross section aL. The solid curves are calculated in the laboratory frame, 
whereas the dotted curves in the cm. frame. The triangles are the real photon points | |26| . 
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FIG. 3. The separated longitudinal and transverse differential cross sections of the reaction 

and 



3, e'TT"*")?! are plotted as function of the virtual photon four- momentum squared k'^ at 



=1175 MeV. The data points are from Saclay [27| and Bonn 
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FIG. 4. Separated cross sections dar/d^Tr + e. da / dO.-,^ as a function of W at fixed A:^ and 6.,^. 
The solid lines are calculated with the full operator; dotted lines with the same form factors and no 
gauge fixing terms; dashed lines with the quadrupole component left out. Data are from Ref. |28|. 
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5. Same as in Fig. H, but for the interference cross section daTh/d^-n 
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FIG. 6. Pion azimuthal angular dependence of the virtual photoproduction cross sections are 
compared with the data [|8| at /c^ = _o.l5 {GeV/cf, W=1160, 1200 MeV and 0^ = 15°, 25°, 35°. 
The legends are the same as in Fig. H. 
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FIG. 7. The 9-,^ distributions for the virtual cross section of p(e, e'Tr^)p are compared to the 
data from Bonn |29| for two pion azimuthal angles. The solid lines are calculated with the full 
operator; dotted lines with the same form factors and no gauge fixing terms; dashed lines without 
the quadrupole component. 
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FIG. 8. Polarization cross section daTr/d^-K 

function of W at /c^ = _i [GeV/cf and 4 
values of 9.,^. The solid lines are calculated with the full operator; dotted lines with the same form 
factors and no gauge fixing terms; dashed lines without the quadrupole component. Data are from 
Ref. pa. 
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FIG. 9. The coordinate system of the reaction A{e, e'7:N)B in the laboratory frame. 
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30 60 90 120 150 180 

(deg) 

FIG. 10. Kinematic variables as a function of the pion angle for one case of Kinematics I. The 
scale for angles is on the left side and the scale for energies is on the right side. The specified 
variables are: u;=400 MeV, e=0.95, k"^ = -Q.lh{GeV/ cf , (p^ = Q, = 180°, Q=100 MeV/c and 
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FIG. 11. Pion angular distributions of the virtual differential cross sections for vr^ electropro- 
duction from ^^C are plotted under kinematics I, along with the contributions from the background 
and the resonant terms in the production operator. The kinematics and the legends to the curves 
are given in the figure. 
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FIG. 12. Sensitivities to the quadrupole component of the A excitation for tt electroproduction 
under kinematics I. 
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FIG. 13. k'^ dependence of the cross sections under kinematics I. Also shown are the sensitivities 
to different choice of form factors. 
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FIG. 14. Distortion effects due to final state interactions of the pion and the nucleon with the 
residual nucleus under kinematics 1. 
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FIG. 15. Pion angular distribution of the differential cross sections under kinematics II are 
plotted along with the contributions from the background and the resonant terms in the production 
process. 
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FIG. 16. Sensitivities to the quadrupole component of the A excitation for tt electroproduction 
under kinematics II. 
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FIG. 17. Pion angular distribution of the differential cross sections under kinematics III are 
plotted along with the contributions from the background and the resonant terms in the production 
process. 



44 




e„ (deg) 





FIG. 18. Sensitivities to the quadrupole component of the A excitation for tt electroproduction 
under kinematics III. 
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FIG. 19. fe^ dependence of the longitudinal and transverse cross sections in the direction of the 
virtual photon is shown under kinematics III along with sensitivities to the form factors. The full 
curve is with realistic form factors, the short-dashed curve is with dipole pion form factor instead 
of monopolc, the long-dashed curve is with all form factors equal to F{k'^) = 1/[1 -|- (-A;2)/0.71]2, 
and the dot-dashed curve is with no form factors. 
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FIG. 20. Out-of-plane distributions of the differential cross sections are sliown as a function of 
the average azimuthal angle (f) = {(f)^^ + 4>p)/2 at fixed = ^tt — = 135° using kinematics III. 
Also shown are the sensitivities to the quadruple moment. 
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FIG. 21. Same as in Fig. 20, but at fixed A(/> 



180°. 
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FIG. 22. The dependence of the calculated cross sections of the reaction CO{e,e'TT~p) 
are compared with the data from SLAG §. The dashed curve is the integrated PWIA result 
carried out according to Eq. (|4^. The solid curve is the DWIA result calculated according to the 
approximation given by Eq. (|4^). 
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